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ABSTRACT 

We present observations made with the 10m Heinrich Hertz Sub-Millimeter Telescope of HCN (3-2) 
emission from a sample of 30 nearby galaxies ranging in infrared luminosity from 10 10 — 10 12 ' 5 L© 
and HCN (3-2) luminosity from 10 6 — 10 9 K kms" 1 pc 2 . We examine the correlation between the 
infrared luminosity and HCN (3-2) luminosity and find that the best fit linear regression has a slope 
(in log-log space) of 0.74 ±0.12. Including recently published data from Gracia-Carpio et al. tightens 
the constraints on the best-fit slope to 0.79 ±0.09. This slope below unity suggests that the HCN (3-2) 
molecular line luminosity is not linearly tracing the amount of dense gas. Our results are consistent 
with predictions from recent theoretical models that find slopes below unity when the line luminosity 
depends upon the average gas density with a power-law index greater than a Kennicutt-Schmidt index 
of 1.5. 

Subject headings: galaxies: evolution and ISM and starburst — ISM: molecules — submillimeter 
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1. INTRODUCTION 

For decades, it has been known that the star for- 
mation rate (SFR) in galaxies is intimately r elated to 
the g as reservoir from which stars are formed ([Schmidt 
1959). Observations of the galactic-averaged surface 
density of HI and CO gas indicate that the SFR in- 
creases with total gas surface density (HI + H2 mea- 
sured by CO J=l-0) according to Esfr ~ E^Lj where 
N = 1.4 ± 0.15 (Ken nicutt-Schmidt law, hereafter KS 
law: lKennicuttlll998bl ). Recent studies have focused on 
dense molecular gas tracers (e.g., HCN (1-0), CO (3-2)) 
and have found a tight, linear correlation between the 
SFR, traced by infrared (IR) luminosity (Lir), and the 
mass of dense gas, traced by m olecular line luminosity , 
L mol (|Gao fc Solomon! l2004aT )hl: iNaravanan et alJl2005h . 
The surface density relation and the SFR-L mo i relation 
(for dense gas tracers) appear to predict different behav- 
iors for the underlying star formation law in galaxies. 

The linear relationship between Lir and HCN (1- 
0) luminosity over th ree decades in Lir found by 
iGao &: Solomon! (|2004al |bl) in local star forming galax- 
ies has been interpreted as a constant star formation 
efficiency (SFR/ M112) traced by dense molecular gas. 
IWu et all (|2005[ ) extended the observed linear relation- 
ship between Lir and HCN (1-0) luminosity to Galactic 
clumps, positing that if HCN (1-0) emission faithfully 
traces dense molecular core mass (above a cutoff lumi- 
nosity of Lboi > 10 4 ' 5 Lq), then constant SFR per unit 
mass is a result of a dense molecular clump comprising 
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a "fundamental unit" of star formation. The observed 
extragalactic linear correlation is a natural extension of 
the constant SFR per unit mass obs erved toward dense 
molecular clumps i n the Milky Way (|Plume et alJl!992l 
IShirlev et al.l l2003. 2007]). In this picture, ultraluminous 
IR galaxies such as Arp 220 which lie on the linear SFR- 
HCN (1-0) relation simply contain more cluster forming 
units, and a higher fraction of de nse molecu l ar gas . 

The interpretation outlined in IWu et all (|2005f) pre- 
dicts a linear relation between SFR and tracers of even 
higher critical (n cr it) density than H CN (1-0). How- 
ever, recent theoretical models from INaravanan et al.l 
(|2007l hereafter N07) and lKrumholz fc Thompson! (|2007l 
hereafter KT07) predict that the power-law index be- 
tween the SFR and higher critical density tracers such 
as HCN (3-2) should in fact be below unity. The physi- 
cal explanation for this behavior is that in systems with 
predominantly low density gas, emission from high crit- 
ical density lines originates in the extreme tails of the 
density distribution resulting in a L mo i-<n> relation 
with a slope greater than unity. This effect drives a 
SFR-L mo i relation with a slope below unity for trac- 
ers with higher n cr i t than that of HCN (J=l-0). To 
test this prediction, we have measured the HCN (3-2) 
line luminosity (LJi CN ( 3 _ 2 )) from a sample of 30 galaxies 
and compared our results with recently published data 
from lGracia-Carpio et ail (|2008t hereafter GC07). Both 
datasets show a Lir - L, 



HCN(3-2; 



slope that is significantly 
below unity, in agreement with the model predictions. 

2. OBSERVATIONS 

Observations of HCN (3-2) {v rest = 265.886431 GHz) 
were obtained from 2007 February through 2007 June 
using the 10m Heinrich Hertz Sub-millimeter Telescope 
(SMT) on Mt. Graham (Arizona). Our sample includes 
30 galaxies and covers a broad range of IR luminosities: 
Lir ~ 10 10 ' 2 — 1O 12 ' 5 L . Central positions of all galaxies 
were observed in HCN (3-2). One nearby galaxy (D < 
7 Mpc) was mapped (NGC0253). The FWHM of the 
SMT at 265 GHz is ~30", such that a single beam covers 
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the central kpc of galaxies beyond 7 Mpc. 

The observations were made using the new ALMA 
sideband separating receiver and the Forbes Fil- 
terbank spectrometer. This dual-polarization re- 
ceiver uses image-separating Superconductor-Insulator- 
Superconductor (SIS) mixers that are significantly more 
sensitive than conventional receiver systems using qua- 
sioptical techniques for image separation. Receiver tem- 
peratures were typically ~ 100 K and system tempera- 
tures ranged from 200 - 400 K. The 2048 channel, 1 MHz 
resolution spectometer was split into two 1 GHz band- 
widths, one for each polarization, corresponding to a ve- 
locity coverage per-polarization of 1130 kms -1 . All of 
our targets were observed in beam-switching mode with a 
chop rate of 2.2 Hz and a beam throw of 2'. We employed 
a 6 GHz intermediate frequency with the HCN (3-2) line 
in the lower side band. Image rejections of the upper side 
band were typically 18-20 dB. 

Our observing strategy involved pointing and calibra- 
tion observations approximately every two hours using 
either Jupiter or Saturn when available, otherwise DR21 
or W30H. Calibration scans were obtained in position- 
switch mode with reference off positions of 5' (30' for 
DR21) in right ascension. We found typical pointing 
errors of 2-3" and measured the main beam efficiency, 
77, n b, to be 0.67 ± 0.04 for filterbank A (H-polarization) 
and 0.80 ± 0.04 for filterbank B (V-polarization). An- 
tenna temperatures are converted to main-beam temper- 
atures using T m b = T^/rjmi,. We assumed a systematic 
calibration uncertainty of ^20%. The integrated inten- 
sity of HCN (3-2) emission was calculated using veloc- 
ity intervals based on the HCN (1-0) line profiles. We 
computed the HCN line luminosity using Equation 1 in 
iGao & Solomon! (|2004b| ) for objects further than 7 Mpc 
away and using their Equation 4 for NGC0253, the 
nearby galaxy which we mapped. Our mapping strategy 
involved 5 pointings along the major axis and 3 along 
the minor axis, each separated by 15". Our results are 
presented in Table [TJ 

3. RESULTS 

In Figure [TJ we show the Lir - L^ CN ^ 3 _ 2 ^ relation us- 
ing the data from Tableland from GC07. Although in 
principle both star formation and Active Galactic Nu- 
clei (AGN) processes heat the dust in a galaxy, the IR 
luminosity from our sample of nearby galaxies is largely 
uncontaminated by AGN, and so we use Lir as a proxy 
for the SFR. Some exceptions are two objects presented 
in this paper (NGC1068 and NG C7469) as well as two 
objec ts from lGracia-Carpio et alj (Mrk231 and Mrk273; 
2008). Excluding these sources from the fit does not 
significantly alter our results (a more detailed investiga- 
tion into the role of AGN will be presented in a com- 
panion paper; S. Juneau et al. in prep). Furthermore, 
these results (particularly our best-fit slope values) do 
not change significantly if one extrapolates from the FIR 
luminosity to estimate Lir. We conservatively assume a 
20% uncertainty in the correction factor needed to gen- 
erate a Ltr value from th e IRAS flux densities (Tab. 1 in 
ISanders fc Mira bel 199(f), which ends up dominating the 
total error budget. Adjusting this uncertainty from 10% 
to 30% does not have a significant effect on the resulting 
best-fit slope. 

We use the publically available Bayesian Monte Carlo 
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NGC0253 


13.67=b2.85 


8±2 


0.25±0.05 


NGC0520 


0.55±0.13 


14±3 


0.8±0.2 


NGC0660 


1.19±0.29 


4±1 


0.31±0.06 


NGC0695 


< 0.13 


< 56 


4.0±0.8 


NGC1068 


6.01±1.22 


38±8 


1.0±0.2 


NGC1614 


< 0.21 


< 23 


3.5±0.7 


NGC2146 


0.68±0.15 


2.6±0.6 


1.1±0.2 


NGC2903 


0.61±0.16 


0.9±0.2 


0.15±0.03 


JNGG3079 


2.58±0.53 


16±3 


0.6±0.1 


NGC3628 


1.04±0.22 


3.6±0.8 


0.18±0.04 


NGC3690 


0.83±0.19 


39±9 


7.1±1.4 


NGC3893 


< 0.23 


< 1 


0.16±0.03 


NGC4414 


< 0.19 


< 0.5 


0.35±0.07 


NGC6240 


0.88±0.19 


225±49 


7.0±1.4 


NGC6701 


0.36±0.10 


27±8 


1.1±0.2 


NGC7331 


< 0.16 


< 0.5 


0.37±0.07 


NGC7469 


0.57±0.14 


65±16 


3.3±0.6 


NGC7771 


0.75±0.17 


66±15 


2.2±0.4 


UGC051017 


< 0.11 


< 74 


9.6±1.9 


VIIZW31 
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a lower limits due to spatial undersampling b mapped 

Markov Chain routines of iKellvl (|2007f ) to compute the 
linear regression between log(LiR) and log(Lg CN ^ 3 _ 2 -j). 
This routine assumes that the distribution of the inde- 
pendent variable can be well described by a mixture of 
Gaussian functions and accounts for heteroscedastic er- 
rors in both Lir and L^ CN ^ 3 _ 2 ^ . The posterior distribu- 
tions of possible slopes and y-intercepts are sampled. We 
define the best fit using the median slope and intercept 
values. The dotted grey line in Figure Q] is the best fit to 
those galaxies observed by the SMT with >2a detections 
(excluding IC342 and M82, for which we have only lower 
limits on L^ CN ^ 3 _ 2 j) and is described by the following 
equation: 

logL IR =(0.74±0.12)xlogL H cN (3 - 2 ) + (5.7±0.9). (1) 

A slope of unity is inconsistent with the data at the 98% 
confidence level. However, we caution that the formal un- 
certainty is large enough to encompass a slope of unity at 
the 2cr level. We note that using a simple least squares fit 
routine (SVDFIT in IDL) instead of the linear regression 
code from Kelly ( 2007) produces a slope of 0.73. 

We can expand our sample by including data from 
GC07. These authors observed HCN (3-2) emission from 
a sample of 13 galaxies with the IRAM 30m telescope, of 
which 10 were significantly detected. We observed 9 of 
these sources and were able to detect 6 of them. Com- 
bining our sample with that from GC07, we can place an 
even stronger constraint on the value of the slope. For 
the sources marked in green triangles, we use the aver- 
age Lhcn(3-2) va l ue 0:l the two surveys (measurements 
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Fig. 1. 



Lir as a function of Li 



HCN(3-2) 



for the sample of 



extra-galactic sources reported here (open circles), those observed 
by GC07 (blue squares), and those detected by both surveys (green 
triangles). Objects not detected above the 2a level are shown as 
upper limits while sources within 7 Mpc that were not mapped 
(IC342 and M82) are shown as lower limits. Error bars represent 
la uncertainties. The solid line is a best fit linear regression to 
the galaxies with detections from the combined sample and has 
a slope of 0.79. The dotted and dashed lines represent the same 
quantities but as applied to the SMT sample only and GC07 sample 
only, respectively. The inset shows the distribution of slope values 
returned by the fitting routine for each of these three samples. In 
the combined sample, a slope of unity is ruled out at the 99% 
confidence level. 

of the overlapping objects were in agreement to within 
20%), but the Lir value as given in GC07, since they use 
a more accurate prescription for computing the total IR 
luminosity. 

We note that a similar sample of nea rby galaxies has 
been observed bv lPaglione et al.1 (|1997l ). Unfortunately, 
a direct comparison between their sample and ours is 
problematic because they employed a dual sideband re- 
ceiver system and did not measure the sideband gain 
ratio, which can vary by a factor of two. We examined 
the integrated intensities of the seven galaxies appear- 
ing in both samples and found the values to be consis- 
tent within this factor of two. HCN (3-2) luminosities 
have been publishe d for the nuclear region of a set of 
nearby galaxies bv iKrips et al.1 (|2008f l. However, these 
observations were conducted with the IRAM 30m tele- 
scope (with a beam-size a factor of 3 smaller than the 
SMT) and therefore provide only lower limits to the full 
L<hcn(3-2) emission. For the purposes of this paper, the 

best comparison sample is that of GC07, so we place our 
focus there for the remainder of the discussion. 

We compute the distribution of best-fitting slopes 
when considering only the published detections in GC07. 
This is shown with the grey dashed line in the inset of 
Figure [U where the median and standard deviation are 
0.63 ± 0.20. The best fit linear regression is shown with 
the dashed line. The next step we take is to combine 
all available data and recompute the best fit slope and 
y-intercept. This results in a narrower distribution of 
slope values with a median below unity: 0.79 ± 0.09. 
The distribution is shown as the solid black line in the 
inset of Figure [IJ and the best fit line is shown in the 
full plot. Using the full, combined dataset, a slope of 
unity is ruled out at the 99% confidence level (the me- 



dian value is 2.5er from unity). Using IDL's SVDFIT, 
we find a slope of 0.79 using the full combined dataset. 
One of the objects in our sample (NGC2146) lies signif- 
icantly off the best fit relation, in the sense of either a 
lower L^ CN ^ 3 _ 2 ) value or a greater Lir. To explore the 

extent to which this affects the resultant best fit slope, 
we remove this object from the sample and re-compute 
the slope, finding a larger slope with a smaller disper- 
sion: 0.84±0.07. This remains signicantly below unity 
and within one standard deviation of the result using the 
full sample. Finally, GC07 find evidence for a change in 
the LiR - L^j^^ relation at Lm ~ 10 11 L Q . Restrict- 
ing our sample to galaxies above this IR luminosity, we 
find a slightly shallower best fit slope of 0.64±0.13. 



4. DISCUSSION 



L' 



The sub-unity slope observed in the Lir - ^hcn(3-2) 
relation presents a challenge to our understanding of 
the molecular star formation rate law. If the constant 
SFR per unit dense gas mass applies to observations of 
HCN (3-2) emission, then a linear relationship between 
Ltr and L{jcn(3-2) s h° u kd be observed. One underlying 
assumption is that the HCN (3-2) molecular line lumi- 
nosity is linearly tracing the dense gas mass. This may 
not be valid when observing unresolved galaxies that in- 
clude large quantities of sub-thermally excited gas. Re- 
cent theoretical models from N07 and KT07 suggest that 
the index between the SFR and higher critical density 
tracers such as HCN (3-2) should in fact be below unity. 
N07 used hydrodynamical simulations of isolated galax- 
ies and equal-mass galaxy mergers coupled with a 3D 
non-local thermodynamic equilibrium radiative transfer 
code to probe the relationship of the dense molecular 
gas phase in galaxies with HCN and CO emission across 
a variety of rotational transitions. Meanwhile, KT07 
use escape probability radiative transfer simulations cou- 
pled with models of turbulence-regulated giant molecular 
clouds. While their work did not specifically explore the 
potential relations between SFR and higher critical den- 
sity lines such as HCN (3-2), fundamentally their con- 
clusion behind the physical driver of the Lir-L^ cn ^ 1q ^ 
relation is similar to that of N07, and thus similar results 
for higher lying transitions are expected. 

According to N07 and KT07, the SFR-HCN (u-l) in- 
dex is parameterized in terms of how the molecular line 
luminosity, L mo i is related to the mean molecular gas 
density, <n>, of a given galaxy. If L mo i grows as <n>, 
the SFR-L mo i relation will have an index of ~1.5 (e.g. 
the observed Lir-Lqq^^ index). On the other hand, 
a L mo i-<n> power-law index equal to the KS index pro- 
duces a linear SFR-L mo i relationship. In this picture, 
the fundamental relationship is the volumetric version of 
the KS relation, and the observed linear SFR-HCN (1- 
0) relationship in galaxies results from the HCN (1-0) 
luminosity on average being related to the mean molec- 
ular gas density by an index similar to the KS index 



(i.e. <n> 1 - 5 ). Since L' 



hcn(i-o) r i ses linearly with 
Lir, HCN (1-0) emission is a useful proxy for the total 
SFR, once the effects of AGN have been properly taken 
into account. 

Alternatively, if the L mo i-<rt> power-law index is 
greater than the KS index, the SFR-L mo i index will be 
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Fig. 2. — Predicted slopes in log(Lus) - l°g(LHCN(3— 2) ) space 
as a function of J transition of HCN, as predicted from theoretical 
models by N07 (grey horizontal lines) and KT07 (yellow shaded 
area at the HCN (3-2) transition). Observational constraints on 
the slope measured from HCN (1-0) and HCN (3-2) emission are 
shown with their error bars. The best-fit slopes are consistent with 
the model predictions from both N07 and KT07. 

below unity. Physically, this can be understood as result- 
ing from high critical density lines originating in extreme 
tails of the density distribution and causing a steepen- 
ing of the L mo i-<n> power law index compared to lower 
critical density lines. N07 predict that the SFR will be 
related to the HCN (3-2) luminosity from galaxies by 
an index of ^0.7. Figure [2] is a reproduction of Fig- 
ure 8 in N07, but with new constraints our data — as well 
as the data from GC07 — place on the SFR-HCN power 
law index. The scatter in the predicted slope is repre- 
sented by the horizontal grey lines and is computed by 
randomly drawing a sample of 19 galaxies (which is the 
size of our combined dataset) out of a set of 100 model 
galaxies 1000 times. The yellow shaded region shows 
the range of expected results from KT07 utilizing their 
publically available code with parameters appropriate for 
"normal" , "intermediate" , and "starburst" galaxies. In 
the KT model, the SFR-HCN (3-2) relation is predicted 
to have a slope below unity in the range 0.77 - 0.93. The 
plotting symbols show the best-fit and la range in slopes 
for t he datasets pres e nted he re as well as that of GC07 
and iGao fc Solomon! (|2004bD . All HCN (3-2) measure- 
ments are consistent with the models shown here. 

It is important to note that we do not account for po- 
tential chemistry-related issues. While some evidence 



may exist for HCN-related chemistry in the vicinity of 
the hard X-ray flux associated with an AGN, both theo- 
retical a nd observational inve s tigations have found mixed 
results dLintott , k Vitil 120061; iGracia-Carpio etHl 120061; 
iMeijerink et al.ll2007l : ICombesll2007D . Studies ofthe nu- 
clear region of galaxies like that done by iKrips et alJ 
(|2007f) ofthe Seyfert 2 galaxy NGC6951 will be essential 
to understanding all of the potentially relevant chemi- 
cal processes. Finally, our sample includes only nearby 
galaxies; the most distant objects in our sample lie 
roughl y 100 Mpc away ( Arpl 93 and NGC6240). A mong 
others. IGao et all (|2007l ) and lRiechers et all (|2007h have 



found evidence that the Lir - L, 



HCN(3-2; 



relation steep- 



ens at high redshift and/or high Ltr . 



5. CONCLUSIONS 



We present observations of HCN (3-2) emission from 
a sample of 30 nearby galaxies ranging in IR luminosity 
from 10 10 - 10 12 - 5 Lq and HCN (3-2) luminosity from 
10 6 — 10 9 K km s _1 pc 2 . We find a best-fit slope in 
log(LiR)-log(L^ CN j 3 _ 2 ^) space of 0.74 and exclude aslope 

of unity at the 98% confidence level (although the formal 
uncertainty is large enough to include a slope of unity at 
the 2a level). Adding data recently published in the liter- 
ature yields a slope of 0.79 and tightens the distribution 
of possible slopes such that a slope of unity is excluded at 
the 99% confidence level for this sample of galaxies. This 
sub-unity slope may be an indication that the HCN (3-2) 
molecular line luminosity is not linearly tracing the dense 
gas. Our results are consistent with predictions from re- 
cent theoretical models by N07 and KT07, who predict 
slopes less than unity when the line luminosity - average 
gas density relation has a power-law index greater than 
the KS index. We wish to emphasize that the results 
shown here are pertinent to HCN (3-2) molecular line 
emission and do not contradict results from previous ef- 
forts showing a tight, linear correlation between Ltr and 
Lihcn(i-o)' Indeed, the models discussed in this paper 
successfully account for this behavior as well. 
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